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We present the theory of interaction between a polar vibration in a semiconductor and an electromagnetic 
mode of a surrounding cavity. Tuning the cavity frequency near the transverse optical phonon frequency couples 
the phonon-induced polarization field within the dielectric and the electromagnetic field in the cavity. Depending 
on engineering parameters, we predict that this cavity quantum electrodynamic interaction may reach the strong 
coupling regime. The resonances of the system are in the terahertz spectral region, and while spectroscopic 
measurements are a possible route for the detection of such a system, we emphasize the possibility of measuring 
the Rabi oscillation directly in the time domain with femtosecond optical pulses. 

PACS numbers: 71.36.+c,42.50.Pq 



I. INTRODUCTION 

Weisbuch et al. initiated a new field of study with 
the demonstration of Rabi splitting between excitons 
and electromagnetic modes in nanoengineered semi- 
conductor heterostructures , and similar phenomena 
have now been demonstrated in quantum-dot systems : . 
Fundamental and technologically promising develop- 
ments are ongoing, including prospects for application 
to light-emitting diodes, probes of the quantum statis- 
tical limit , and investigation of bosonic many-body 
effects 4 . 

Coincidentally, time-resolved studies of fast inter- 
actions in semiconductors have emerged with fem- 
tosecond laser technology 5 ' 6,7 . With new technol- 
ogy allowing time-domain, optical observations of 
fundamental excitations in crystals, coherent optical 
phonon oscillations have been read out in several sys- 
tems, such as GaAs and bismuth . Experiments have 
reported intriguing light-lattice interactions, such as 
phonon squeezing 10 , pulse-controlled spatiotemporal 
lattice coherence 1 , and exotic phonon dynamics 12 . 

We are proposing a measurement which may reveal 
a Rabi splitting between a coherent optical phonon in 
a polar crystal and a far-infrared cavity mode. We 
mean "coherent" phonon to be taken in the contempo- 
rary sense, referring to a mode exhibiting time-domain 
phase and amplitude information. If such a Rabi split- 
ting can be probed through a time-domain quantum 
beat, the phenomenon may be interpreted as a co- 
herent oscillation of energy between a polar mechan- 
ical mode in the crystal and a vacuum cavity reso- 
nance. As such, the system we suggest may demon- 



strate behavior which is interesting from the standpoint 
of narrow-band terahertz technology and quantum op- 
tics, and which presents favorable engineering param- 
eters. For instance, as will be discussed, a cavity-GaAs 
apparatus should exhibit a Rabi splitting of approxi- 
mately 30 K. This value roughly corresponds to the 
bulk longitudinal phonon frequency enhancement and 
the Kramers-Kronig-accompanied transverse phonon 
oscillator strength. This robust Rabi coupling may 
be accessible with the relatively liberal requirement of 
several-micron-scale engineering. 

The Rabi interaction we are postulating is formally 
and physically similar to time-domain beating ob- 
served in two other coherent phonon systems. At 
high pump fluences, a plasmon-mediated beating has 
been demonstrated , and more recently the phonon 
has been Rabi-coupled to superlattice-tuned excitonic 
states . The intrinsic lifetime of the phonon mode and 
our estimate for the lifetime of the cavity mode are each 
greater than the fast and fluence-dependent dephasing 
time of the plasmon, and than the excitonic transition 
lifetime. We shall discuss our view that the phonon and 
cavity-mode lifetimes will limit the coherent oscilla- 
tion we propose. For the photon mediated case we are 
addressing, this suggests a more pronounced demon- 
stration of strong-coupling behavior, and a greater du- 
ration of the oscillation. 

With the above in mind, the qualitative features of 
our model system take shape. The model possesses 
an in-plane translational invariance, so that the pump- 
driven optical phonon should maintain phase over in- 
plane lengths at least comparable to the dielectric thick- 
ness. The effectively one-dimensional system remain- 
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ing consists of a nodeless electromagnetic cavity mode 
superposed upon a uniform polarization, which fol- 
lows a step function over the dielectric length. The 
polarization is modeled as the combined effect of N 
charged harmonic oscillators, one for each GaAs ionic 
pair. Each oscillator has a natural frequency and a 
mass corresponding to the GaAs reduced mass, and a 
charge corresponding to the Born-effective charge of 
bulk GaAs. 

Phononic oscillator strengths have been phenomeno- 
logically understood for many years 15 ' 16 , and it is no 
surprise that a coherently excited phonon may radi- 
ate. Terahertz radiation from pump-excited coherent 
phonons has indeed been observed 1 . The radiation 
spectra include an emission peak at the longitudinal- 
optical phonon frequency and a broadband response 
from the transient dynamics of photocarriers 1 °, We are 
suggesting the possibility that a coherent phonon, in 
proximity to a cavity of well-defined resonance, can 
undergo a hybridization of its quantum state with the 
electromagnetic resonance, as is the case in atomic and 
excitonic systems. In the weak-coupling limit, sponta- 
neous photon emission and absorption are influenced, 
and in the strong-coupling limit emission and absorp- 
tion can occur in coherent succession. It may then 
be possible to control and enhance the amplitude and 
phase coherence of the excitations. In fact, it has now 
been experimentally confirmed that engineered geome- 
tries can enhance the terahertz radiation of coherent 
phonons by two orders of magnitude 1 and that the 
polaritonic spectrum can be tuned by varying electro- 
magnetic boundary conditions . 

If the phenomena we describe can be observed, and 
optimally controlled, they may provide a promising 
specimen for the investigation of macroscopic quantum 
phenomena. They may also prove relevant to important 
new terahertz technologies. For instance, the system 
we describe may allow more detailed study the phonon 
signal limiting electro-optic sampling, and is related 
to the control of terahertz polaritonic waveforms and 
phonon-polariton photonics- 1 '—. The amplitude decay 
may also help to illuminate the nature of the coherent 
state, and probe interactions among lattice, electronic 
and electromagnetic degrees of freedom. 

We will refer to the coupled phonon-vacuum reso- 
nance excitation as a phononic cavity polariton. Tra- 
ditionally, the polariton was considered as a bulk 
excitation , but there is a conceptual resemblance be- 
tween the canonical polariton dispersion and a cavity- 
resonant system. If one imagines a level (avoided) 
crossing diagram representing the variation of the me- 
chanical and electromagnetic energies with inverse 
cavity dimension, the result is not unlike the standard 
polariton dispersion, with the correspondence of wave 
vector to inverse cavity dimension. In this picture, 
the mechanical resonance is constant at the transverse- 
optical phonon frequency, while the electromagnetic 



resonance energy increases linearly with the inverse di- 
mension. Thus Fig. 2 bears resemblance to the polari- 
ton dispersion of bulk GaAs. 

Accordingly, Savona et al, in a study of excitonic 
cavity polaritons, generalized Hopfield's bulk polariton 
theory to allow for a vacuum cavity, and for the exten- 
sion of the electromagnetic field beyond the boundaries 
of a confined condensed-matter system- 3 . In this work, 
we develop the analogous theory for the phononic cav- 
ity polariton. 

II. DERIVATION OF INTERACTION 
HAMILTONIAN 

The formalism to follow will emulate Savona's light- 
matter interaction theory for exciton polaritons, which 
motivated it. Our system is comprised of a rectangular, 
semiconducting crystal, of thickness d c i, inside a cav- 
ity bound by far-infrared reflective surfaces. The half- 
wave resonance length, including the dielectric and the 
cavity, is d c . In realizing such a system, a convenience 
is presented by the high reflectivity of simple metals in 
the far infrared. This property makes the use of Bragg 
reflectors and other advanced micro-fabrication tech- 
niques, which are required at optical frequencies, un- 
necessary in the far-infrared. The semiconductor pro- 
viding lattice motion must lack an inversion symmetry, 
so that it possesses an optically active phonon. The 
simplest, and best studied such material is GaAs, and 
so we will adopt it in the theory to follow. The geome- 
try and relevant lengths are illustrated in Fig. 1. 

The far-infrared frequencies are much smaller than 
the fundamental electronic gap of GaAs, and so, ne- 
glecting any photocarriers, our Hamiltonian need only 
consider the bare piece Hq which includes the phonon 
field in the lattice and the electro-dynamic field, and 
their interaction, i.e. 

H=H +H l . (1) 

Generally, the lattice term is a set of harmonic oscilla- 
tors, each with a wavevector within the Brillouin zone, 
and there are three branch index values per atom in the 
primitive unit cell. Similarly, the photon field can be 
expanded into components of given wavevector and po- 
larization, such that the reflective boundary condition is 
maintained. In the crudest approximation, we limit our 
consideration of the electro-dynamic degrees of free- 
dom to the lowest-energy mode satisfying the reflective 
boundary condition, the half-wave resonance. The gen- 
eral multimode case is treated in the Appendix. Fur- 
thermore, we limit the lattice-light interaction to lat- 
tice modes of vanishing crystal momentum, relative to 
Brillouin zone dimensions. This is justified from en- 
ergy and momentum conservation in the emission and 
absorption of photons by lattice vibrations 24 . Finally, 
both the phonon and photon polarizations are taken to 
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FIG. 1: Schematic of a dielectric crystal with polarization 
current P = Py inside a cavity. The length scales satisfy d c ~ 
cn/2(i>jo and dj < d c <C l x ,l y where l x , Y are the transverse 
dimensions of the dielectric. 



be in the i-direction, which is within the plane of the 
interface and reflective surfaces. 

With these approximations in mind, the photon and 
phonon contributions can be written explicitly in terms 
of the field creation and annihilation operators, 



H = rm TO (a J( a+ 1/2) 



Tlftc 
F 



b + 1/2) + H h (2) 



The first term on the right hand side is the lattice term, 
and indicates the transverse optical phonon frequency. 
The second term is the photon term, and indicates the 
speed of light and the effective dielectric constant e, 
which characterizes the composite system of cavity and 
dielectric. A reasonable approximation to e, which will 
be assumed here, is a volume average of the bulk semi- 
conductor dielectric constant and the permittivity of 
free space; i.e. 



£ = ! + (£()- lK,/£/,. 



(3) 



A more accurate value may be obtained through meet- 
ing Maxwell's boundary conditions at the interfaces 
and reflective surfaces, but such a correction only shifts 
the cavity resonance, and we are tuning near the trans- 
verse optical phonon frequency, COro, in any event. In 
the second term on the right, the creation and destruc- 
tion operators correspond to photons with polarization, 
again, along the i-axis. 

The derivation of our interaction Hamiltonian 
closely resembles Savona etal's excitonic theory. The 
general form of the interaction, with electromagnetic 
vector potential A polarized along x, is: 

1 72 



-z 



Hi = —VA(r)-[x(r),H p ] 



ih 



2/jc- 



This sum is taken over the lattice site positions, where 
each Ga and As atomic pair are interpreted as har- 
monically oscillating dipoles, characterized by Born- 
effective charge, Z, reduced mass, /j, and frequency 



(Oto- The x-projection of the normalized, local, har- 
monic oscillator displacement field can be expressed 
through the phonon operators above: 




{a + a f ) 



(5) 



where there are N atomic pairs within the dielectric. 
The commutator in the interaction Hamiltonian, then, 
is: 



x(r),H p 



' ft&TO 



(a - 



(6) 



Now, the normalized vector potential is approxi- 
mated as a planar, standing halfwave, polarized along 
x. The phase gradient is along the interface normal, 
z, and is approximated by its length-averaged value, 
q = n/d c . If the dielectric and the cavity share a cross- 
sectional area S, then, the i-polarized component of the 
vector potential is: 



A(r z ) = 



2%hc y /2 r 



(i) 



The first term in the interaction Hamiltonian may be 
re-written as: 



■■■ZrA{r)^[x{r),H v ] = fl r (g|) 



1/2 



fift)T 



2N(i 

cos[gr,](a — a' i ')(b + Z? t ). 

(8) 

The sum over atomic pairs, each occupying a volume 
£1, may be approximated with an integral over the z 
direction: 



S f d d/ 2 

VcosIoTt] = — / dr?cos[qr-] 

y " &J-d d /2 



(9) 



The final form of the first term in the interaction Hamil- 
tonian follows, 

(10) 

Where NO. = Sd c i, and the conventional definition of 
the nuclear plasma frequency is substituted, C0 ; , = 
Z(47i/ ( ui2) 1 / 2 . The second term in the interaction 
Hamiltonian is: 



L r A 2 (r) = L- 



2nhc 



2/jc 1 Sqd c </£ 



■[qr^{b + tff 



dc 



(b+b^) 2 . 

(ID 

With these explicit expressions for the interaction 
term in the Hamiltonian, the whole Hamiltonian for the 

system is, 

H = h& T0 (a j 'a+ 1/2)- 



+ih(i>Ap(b + b 1 ) (a — a ' ) 
+h(aAA(b + b r ) 2 . 



(12) 
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The interaction coefficients are taken from above, 




Physically, this interaction corresponds to a uniform 
polarization current coupling to the positive and nega- 
tive wavevectors of the lowest frequency standing wave 
that matches the boundary conditions. 

IH. NEAR RESONANT BEHAVIOR 

Simplistic physical considerations lead to the con- 
clusion that any eigenstate of our system requires a 
mixture of lattice and electromagnetic character. A 
phonon will radiate, and conversely a photon will be 
absorbed by the crystal, so that neither one of these 
states alone is stationary. To find an eigenstate of the 
Hamiltonian, which we refer to as the cavity polari- 
ton, we define its operator subject to two conditions. 



First, we require the polariton to be some superposition 
of phonon and photon fields; and second, we require 
that the commutator of the polariton operator with the 
Hamiltonian take the standard harmonic form. 

The former is satisfied by writing the polariton field 
operator as a linear combination of creation and anni- 
hilation operators of the phonon and photon fields, 

B = Wa+Xb + Ya f + Zb f . (14) 

Each of the coefficients, W, X, Y, and Z, implicitly de- 
pends on the cavity dimension, d c , through the cou- 
plings, Ti(£>ap and /iCOaa- The second requirement is ex- 
pressed as, 

[B,H]=EB, (15) 

where on the right-hand side the energy eigenvalue is 
indicated. 

Making the substitutions, D = Hg>aa, and C = h(£>Ap, 
the above two equations may be combined in matrix 
form: 



hue 



2D-E 



iC 
2D 

-iC 



-iC 

h&TO — E 
-iC 





line 



-2D 
-iC 

-2D-E 



iC 



-iC 


-iC 
-h&TO — E 



/ W 
X 
Y 



The eigenvalues are obtained by solving for the condi- 
tion that the above 4x4 matrix have zero determinant. 
The exact expressions for positive branch of energies 
that exhibit an anticrossing is: 

E± = h (co 2 +4coaaK) + cOj ±(co 4 + 8coaaC0 3 
+ 16(0^4 co 2 - 2(0 2 - O o> 2 - 8co^O)j a) 

+ 16(0^0)7-0(0 + 03 4 TO ) l/2 )/ 2 

(16) 

where (0 = cn/y/£,d c - These energies, as a function of 
d c , are plotted in Fig. 2. When the cavity is tuned far 
from (Oro, the two branches have negligible curvature, 
with the slope of the photon-like branch proportional to 
c, and the slope of phonon-like branch approximately 
zero. In this limit, the cavity mode and lattice vibration 
are decoupled. When the cavity mode is tuned near 
(Oro, however, there is a significant coupling, and there 
is evident curvature in the two branches. 

The Rabi frequency at maximum coupling can be 
seen in Fig. 2 as the minimum value of the energy 
difference between the two branches. The figure also 
indicates the cavity detuning range over which signif- 
icant coupling may be expected. The parameters for 



GaAs are (£>to = 36 microns, (0 ; , = 51 microns, and 
£o = 10.89. The crystal width to cavity width ratio is 
taken to be fixed at d c \ = d c /5Q. 



IV. DEPHASING AND PROSPECTS FOR 
MEASURING RABI OSCILLATION 

When a cavity-mode is weakly coupled to an excita- 
tion of non-zero oscillator strength, spontaneous emis- 
sion and absorption rates can be influenced by tun- 
ing the cavity through the resonance. When the cou- 
pling is increased, and the inverse energy difference 
of the avoided crossing surpasses a characteristic dis- 
sipation rate, coherent interaction may be observed. 
For instance, a Rabi oscillation can take place between 
modes of primarily lattice character, and of primar- 
ily electromagnetic character. The fundamental dis- 
tinction between the strong- and weak-coupling con- 
ditions is governed by the relative magnitudes of the 
near-resonant oscillator strength of the excitation, and 
the mode and excitation lifetimes. 

Near cavity resonance, the strong coupling condition 
is coap > |Yc — Yro|/4, where JcJto are the full-width 
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FIG. 2: Variation of the positive energy eigenstates of the po- 
lariton Hamiltonian as a function of the inverse cavity length 
d c . The dielectric is taken to be a GaAs crystal of width 
4/ = del '50. The resonance occurs near d c ~ ncj '^/ecoro = 
16.45^111. 



at half-maximum linewidths for the cavity mode and 
phonon . For d c \ <C d c this implies 



Wp^dd/d c > \y c -jTo\/2. 



(17) 



The bulk GaAs phonon linewidth is well-studied" 1 , 
and is in the thousands of microns. The cavity lifetime 
is a less fundamental quantity, and depends on engi- 
neering of the cavity. State-of-the-art etalon techniques 
can increase the cavity lifetime at the expense of the 
tuning range. With a two-micron tuning range, a cavity 
lifetime which is again in the thousands of microns is 
possible with commercially available technology. We 
see than that the above inequality is satisfied to an order 
of magnitude. This strong coupling is a consequence of 
the two-dimensional emitter at the interface. 

We conclude then that strong coupling is likely. The 
Rabi period is also greater than the temporal resolu- 
tion of ultrafast techniques using femtosecond optical 
pulses. Presuming that the dependence of the reflectiv- 
ity on lattice strain is distinguishable from the effect of 
the electromagnetic field, the Rabi oscillation should 
be observable by pumping the optical phonon, and per- 
forming time-domain reflectivity measurements. 



V. SUMMARY AND CONCLUSIONS 

We express the quantum interaction between an op- 
tically active phonon and a surrounding cavity's res- 
onant electromagnetic mode. When the two are tuned, 
we demonstrate a Rabi splitting, which is evident in the 
calculated energy levels for varying cavity dimension. 
Furthermore, we review various experimental studies 
and analyze the strong-coupling criteria to support our 
view that a time-domain Rabi oscillation may be ob- 
served in the system we describe. 



APPENDIX A 

In the following, we present a full multimode cal- 
culation of the the polariton Hamiltonian and the dis- 
persion relation for a dielectric crystal embedded in a 
cavity. The bare Hamiltonian for our system is (up to 
an overall constant): 



(Al) 



Here a^ a ,al are the annihilation and creation op- 
erators for the phonon field at the frequency &to- 
The phonon wave vector q = (q±,q z ) is taken to be 
small (so that wave lengths are macroscopic), and the 
phonon polarization vector is e a . Similarly, b^ a ,bt g 
are the annihilation and creation operators for the pho- 
ton field. The dielectric constant is the volume av- 
eraged value e = 1 + (eo — \)dd/d c . In the radiation 
gauge, the Hamiltonian describing the interaction be- 
tween the fields is: 



Hi 



H 



d'xA-P- 

crystal 
1) uP) 



2Q./JC 2 J crystal 



d'xA-A 



H 



(A2) 

where £2 is the volume per lattice cell of atomic pairs, 
/j is the reduced mass of the pair, Z is the Born effec- 
tive charge, and c is the vacuum speed of light. The 
macroscopic polarization current is 



r,z 2 



i(q y y+q z Z) 



(A3) 



The vector potential must satisfy appropriate bound- 
ary conditions along the cavity axis £. We 
assume perfectly reflecting mirrors at the cavity 
boundaries meaning the momentum components of 
the field along z satisfy q- = nn/d c for n € Z. 
The suitably normalized vector potential is then 
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A = 



1 2nhc 
^/iSd, 



sin(2«7iz/ii f ) 



z -(| ? l| 2 + (2«7rM.) 2 ) 1/4 

cos((2n+l)7iz/t/ c ) 



(| 9 l|2 + ((2n+l)7 t K)2)l/4 / 



(A4) 



Hereafter we suppress the polarization index o as it is 
understood that only x polarized fields contribute to the 

interaction. Using the relation P = — ^[P,//o], we ob- 
tain the full expression for the interaction terms: 



cos((2n + l)jiz/c/ e ) 



H ' l) = lTmp ^^LJ\^\ (,? + wn + wwf* 



9yi — 3j — * "7vi — — — tyv. 



sin(2imz/d c ) 



(qj+(2nn/d c ) 

EE( 

q y ,q'.neZ 



-'(ly+l'y)x e -iq'-z^j 



C 



2mt fl_, 



o'. +ft„ ImtflJ, . + b t lim a„ xlL +b 



,t „t 



t .i,t „ ._«,t 



(2) 



flag / 3 / cos((2m + \)%z/d c ) cos((2n + \)%z/d c ) 

^^Jcrystal * £ J^($ + ((2m + d c ffl\qj + ((2n + 1)^/^)2)1/4 



x (ft (aB +i), e'fr y + fr' (2 ,„ +1) , e '^)(ft (2,,+D.e 



x (ft„ w^ + ft* 2™^)(ft„ v ^^ + ft; ^e-^) 



'lyy _i_ ht p~^yy\ A sin(2m^/rfc)sin(2n7ic/flf c ) 

1y^ 2 't ' (^+(2'W^) 2 ) 1/4 (??+(2™K) 2 ) 1/4 



„ v 2™ +ft av 2 ™tft"'' 2ra +ft„ 2mnft _, 2™ 2MI ^ „ Inn ) 



E E [p\ qy ^)^ b qy,^ b 

D ( ?y ,(^),^(\ > ^ fo -^P^ + \,,(™%,,<^ ^(2^%,^ 

, . . . ! (A5) r 



)uft' ra-j-iw+ft ft-nk* f^+Dn+ft «_,n.ft 



In the integration we have assumed that £ y 3> kc/(0to Note that the interaction is invariant under the parity 
so that the Hamiltonian is approximately translation- transformation q,q' -»■ -q,-q'. The coupling coeffi- 
ally invariant along the y direction. The momentum cients are related by 



dependent coefficients in Hj corresponding to cou 
pling with even and odd parity photonic fields are 

C e - Vim / m ro ! Y 

^(q y ,qzU'z ~ nW P\/ cd c d ds /l (ql- q >2)(qj+q2y/4 X 

(q z cos(d c /qi/2) sin(d^q z /2) 
-q' z cos(d d q z /2) sia(d d q' z /2)), 



C° 



In the bulk limit, d c i — > d c , and for an infinite dielectric 
the momentum components satisfy \q z \ = \qi\. Then 
we recover the appropriate infinite dielectric coupling 
relations (see e.g. ): 



(q y ,q z ).,q' z ~ nw P\J cd c d d </e ( 9 2- 9 «)(,2+^)l/4 ^ 

(q'-COs(d c iqi/2) sin(d c iq : /2) 
-{q z cos(d d q z /2)sm(d d q' z /2)). 

r(2) 



/ico, 



(A6) 



D 



l6c^/eQ^qj+q% 



(qy,qz),q' : Sc^eSy/qJ+qJ' 



Similarly, for H, the momentum dependent coupling .,„ ... c 

— ■ , ' , , , rL . The total Hamiltonian tor our system is then 

coefficients between even and odd held components is 



D 



ha 1 ,, 



{qy,qz),q'z 4^«/ c {qJ+qfF*(q$+qJ^ 

■ )d d /2) ± sm((q z +q' z )d d /2) 



(A7) 



H = E L h ^\j4y + ( nn / d c) 2b qy , m /d c b qy,™/dc 
q y neZ V c 

1z 

(A9) 
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Since the Hamiltonian is a quadratic polynomial in 
field operators, H can be diagonalized by a Tyablikov- 
Bogoluibov transformation • . The Hamiltonian is 
translationally invariant along y so different q y momen- 
tum components do not mix. To obtain the solution, we 
introduce a new set of bosonic operators 

+E i W Mz)%^ + Y (ly > Jz) a Uy-d ■ 
<?'- 

(A10) 

The coefficients of the expansion are determined by 
satisfying the commutation relation 

[B qy ,H]=EB qy . (All) 

The condition that a solution exists then reduces to the 
following dispersion relation which is an implicit func- 
tion of energy: 

E 2 - {h<» TO f = £ (fh + dE)+f2i(E)), (A12) 
where 

< ) . ) ,+^r L[C^ ) ] 2 ) 

E 2 (hl» p ) 2 c l d /2d c . sm(£nd d /d c ) \ 

e(E ,2 -(fic) 2 (q 2 .+(H/dc) 2 )/z)" indjdc )■ 

(A13) 



In the bulk limit, for each momentum wavevector q we 
obtain the appropriate dispersion relation 16 , 



(*-<*»> 1 )(*-» 




(A14) 



When the cavity is tuned near the resonance condi- 
tion at minimal separation, d c « Kc/(Oto, then the even 
parity q z = ±n/d c modes of the field are strongly cou- 
pled to the polarization current. Restricting attention 
to the q\ = modes of the polarization current which 
dominate the coupling and defining the creation opera- 
tor for the even parity photon mode: 



frt= °'-^V W , (A15) 



and a ' = al , the interacting piece of the Hamiltonian 
then reduces to Eq. 12. 
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